Nonlinear entanglement witnesses 
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Entanglement detection typically relies on linear inequalities for mean vàlues of certain observables 
(entanglement witnesses), where violation indicates entanglement. We provide a general method to 
improve any of these inequalities for bipartite systems via nonlinear expressions. The nonlinearities 
are of different orders and can be directly measured in experiments, often without any extra effort. 
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Entanglement detection is one of the fundamental 
problems in quantum information science. On the one 
hand, it is crucial for experiments, since the question 
whcther a produced state is entangled or not may decide 
whcther a given experiment has been successful or not. 
On the other hand, it is also a challenging task for theo- 
reticians, since the separability problem is, despite a lot 
of progress in the last years, essentially not solved. 

In most of the experiments, so called entanglement wit- 
nesses are used for entanglement verification 0- Q- S d 
0- El • These are linear inequalities for mean vàlues of 
certain observables, if the inequality is violated, the state 
must be entangled. For instance, Bell inequalities are 
such linear inequalities and may hence be viewed as en- 
tanglement witnesses. The question whether one can use 
nonlinear inequalities for entanglement detection has also 
been discussed and many nonlinear criteria are known 

iSEIilElEi. 

In this situation it is natural to ask whether it is pos- 
sible to improve a given linear witness by some nonlinear 
expression. For discrete systems two examples for such 
an improvement are known. Uffink showed that certain 
Bell inequalities can be improved by nonlinear expres- 
sions . Recently, Hofmann and Takeuchi proposed sep- 
arability conditions called "local uncertainty relations" 
[ïfl | , which can improve witnesses in some cases jjjj . For 
continuous variable systems it is known how to summa- 
rize certain famílies of linear inequalities to a single non- 
linear one [ï^ |. 

In this Letter we provide a method which allows to 
improve all entanglement witnesses for discrete bipar- 
tite systems by nonlinear expressions. Surprisingly, this 
works also for optimal entanglement witnesses, that is, 
witnesses where no other linear witness is stronger. Due 
to the simplicity of our method the nonlinear witnesses 
may be used from the beginning, without considering lin- 
ear witnesses anymore. Our method allows to calculate 
such nonlinear witnesses to an arbitrary order. We dis- 



cuss in detail the strength of our constructions for the 
case of two qubits. Interestingly, when implemented in 
an experiment, the nonlinear expressions in our construc- 
tions often require measurement of the same observables 
as the original witness. This allows for an improved en- 
tanglement detection without any extra effort. 

Without loosing generality, we consider witnesses W as 
observables with a positive mean value on all separable 
states, thus a negative expectation value signals the pres- 
ence of entanglement. Our aim is to find a nonlinear wit- 
ness J- as a functional of the type T(g) = Tr(Wg) — X(g) 
which still should be positive on all separable states. 
Typically, the nonlinearity X(g) will be a quadratic poly- 
nomial of certain expectation vàlues. We consider only 
T(q) which are strictly stronger than the witness W. 
That is, we require that T(g) deteets all the states that 
are detected by W and some states in addition. 

Let us explain the main idea for the case of witnesses 
coming from the separability criterion of the positivity of 
the partial transpose (PPT) [lj]. By definition, a quan- 
tum state g shared between Alice and Bob is separable if 
it can be written as a mixture of product states, that is 
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where pi > and X^P» = 1- Tiren, given a quan- 
tum state g = J^ij ki í?*j>*ïK}01 ® I^X'I its partial trans- 
pose with respect to Bob's system is defined as g Ts — 
Y^ij ki Qíj,ik\i}(j\®\k}(l\- If g is a separable state, it can be 
easily seen that the partial transpose is positive, g Tn > 0. 
Thus, if for a state the partial transpose is not positive (g 
is NPTlthe state must be entangled. Indeed, it has been 
shown |j| that for 2x2 and 2x3 systems a state is PPT 
if and only if it is separable, while for other dimensions 
there are also PPT entangled states. 

For any NPT state go we find that g^ B has a negative 
eigenvalue A_ and a corresponding eigenvector \<j)). An 
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entanglement witness for this state is then 



w = \4>){^\ 



Te: 



(2) 



Indeed, due to the identity Tr(XY TB ) = Tr(X TB Y) for 
arbitrary matrices X, Y, we have Tr(g W) = A_ < 
while for separable (and hence PPT) states Tr(gW) = 
{<t>\Q TB \<t>) > holds. Note that the witness in Eq. © is 
by no means specific: since the PPT criterion is necessary 
and sufficient for low dimensions, witnesses of this type 
suffice to detect all states in these systems. Furthermore, 
such witnesses can be shown to be optimal || . 

To improve the witness from Eq. with nonlincar 
terms, first note that a functional like (X)(X') is convex 
(see Lemma 1 in Ref . ) . This implies that a functional 
like Q = (A) — Yli a i{Xi){X}) is concave in the state. 
That is, if g = J2kPkQk is a convex combination of some 
states, then Ç(g) > J2kPkG(Qk)- 

Let us assume that we have chosen A = \<j>){4>\,Xi — 
\(j)}(tjji\ for an arbitrary 1^) and take a separable state 
g. Then g TB is again separable and can be written 
as a convex combination of product states, g TB = 
J2kPk\ a kbk)(akbk\- For \x) = \a k b k ) we have 

Q{\X){X\) = (XIÏMX) ■ [l-^aiixmmx)]- (3) 



=:F(X) 

Thus, if the polynomial P(x) is positive on all product 
states, by concavity the functional Q is positive on all 
separable g Ta . Then, with the chosen Xi, 
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is a nonlinear improvement of the witness W = \(f))((j)\ TB . 
It is important to note that a term like (X Tb )((X Tb )^) in 
Eq. Q is easily accessible in experiments, even if X Tb is 
non-Hermitean. Namely, we can write X Tb = H + i ■ A as 
a sum of its Hermitean and anti-Hcrmitean part, where 
H and A are Hermitean. Then (X Tb )((X Tb )ï) = {H) 2 + 
(A) 2 holds. With this method, we have: 

Observation 1. (a) Let W — \cj))(<fi\ TB be an entan- 
glement witness. We define Xi = \<fi)(ip\ for an arbitrary 
and s(tp) as the square of the largest Schmidt coeffi- 
cient of \ip). Then 



(x TB )((x TB y) (5) 



is a nonlinear improvement of VV. 

(b) If we define Xi = \<fi)(ipi\, i = 1,...,K with an or- 
thonormal basis \ipí), then 
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T^{g) = (\cb)(cj)\ TB ) -^(X^)((Xf B )t) (6) 



Proof. (a) The squared overlap between a state \ip) 
and a product state is bounded by the maximal squared 
Schmidt coefficient 2], thus P( X ) in Eq. © is positive. 
(b) For the we have in Eq. © EiOcIV'iX^lx) = 
Tr(\x)(x\) = 1 which proves the claim. □ 

This Observation provides a whole class of nonlinear 
improvements, one may pick an arbitrary and com- 
pute the corresponding nonlinearity. In special situations 
one can also adjust the improvement: if an experiment 
produces the state g(p) = pg e + (1 — p)g n , which is a 
mixture of an entangled g e and some noise g n , a given 
witness may detect the states only for p > p cr it, i.e., 
Tr(Wg(p CT it)) = 0. Now one can choose the \ip) in a way 
that the nonlinear terms are large at £>(p C rit)j which leads 
to a significant improvement of the witness for the states 
of interest. Concerning the strength of the nonlinear im- 
provements we can state: 

Observation 2. (a) Let W = \4>){(j)\ TB be a witness. 
A state g can be detected by a witness of the type 
from Eq. JSJ if and only if 



< 
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(b) A state can be detected by a witness of the type T^ 2 
from Eq. JBJ if and only if 



< (4>\{q Tb ) 2 ' 



(8) 



holds. In this case, the state is detected by all nonlinear 
witnesses of the type T 1 ^ . 

(c) If Eq. JSJ) is fulfilled, then also Eq. J7J) holds, thus the 
witnesses of the type T^~> are stronger. Furthermore, 
Eqs. (00 are never fulfilled for PPT states. 

Proof. The proof is given in the Appendix. □ 
Let us add that with the same idea also witnesses 
with higher nonlinearities, e.g. of fourth order, can 
be derived. To do so, note first that functionals like 
g = (A)-a(B)(BÏ)-l3(C} 2 (CÏ) 2 are also concave in the 
state. If we choose A = \4>){<f>\,B = = |0}(02|, 

we arrive at a similar expression as in Eq. ©. Now, P(x) 
is a polynomial of fourth order that has to be positive on 
all product states [ï^ . 

To give an example which we will investigate later, let 
us assume that l^i), \ip2) and \<j>) form an orthonormal 
set, where the biggest squared Schmidt coefficient from 
| -01 ) , |^2) is 1/2 and the one from \<j>) is > 1/2. Then, 
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(9) 



is also a nonlinear witness. 



is a nonlinear witness of fourth order. The positivity of 
the polynomial P(x) can directly be checked. 

It should, however, be noted that higher orders do not 
necessarily provide witnesses, which are much stronger 
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than the witnesses with quadratic nonlinearity. The rea- 
son is that we are not trying to approximate the curva- 
ture of the convex set of separable states (in which case 
higher polynomials would be clearly an advantage): in 
our approach, we also require the concave functional to 
be positive on the convex set of separable states which 
makcs the situation more complicated. 

Let us discuss an application to two-qubit systems. A 
genèric optimal witness for a two-qubit system is 

W(a) = | Tb , (10) 

\(j)(a)) = cos(a)|00) + sin(a)|ll). (11) 

It is useful to express all quantities directly in expectation 
vàlues for observables. We use the abbreviations s — 
sin(a), c = cos(a), s 2 = sm(2ce), C2 = cos(2a) and x\x 2 = 
(a x ® <J x ),Ui — (p y <8> 1) etc. here. The witness is then 

(W) = 2( 1 + z i z 2 + S2{xix 2 + y 1 y 2 ) + C2(z 1 + z 2 )). (12) 

To construct improvements, we choose the four vectors 
|Vi/ 2 } = (100) ± JH))/V2, IV3/4) = (|01> ± |10))/V2, 
where the upper signs hold for i = 1,3, and compute the 
corresponding Xi = |0(a))(0i|, which leads to the four 
nonlinear correction terms Xi :— (Xj B )({Xj B )^) with 

^1/2 = ^[(( c± s )i l ± xiX2 ± y^ 2 + ziZ2 > 

+(c T s)(z 1 + z 2 )) 2 + ((c T s)(x 1 y 2 ~ yix 2 )) 2 ] , 

1 2 
<%3/4 = ^ [((c ± s)(x 1 ±i 2 )+(cí s)(z 1 x 2 ± XxZ 2 )) 

+ ((c T s)(yx T y 2 ) + (c ± s)( yi z 2 T Z W2)) 2 ] ■ (13) 

With these nonlinearities, the nonlinear witnesses 
^"W,^ 2 ), and according to Eqs. (011111) can di- 
rectly be written down. Note that, especially for c = ±s, 
several quadratic terms in the Xi can be measured with 
the same measurements as in Eq. (11211 . Thus their im- 
plemcntation requires no extra effort. 

To investigate the power of these new criteria, we con- 
sider the family of states g(a) = (1 — \4>(a)} (4>(a)\ TB )/3, 
which are separable, but they lie on the boundary of 
the set of separable states, since Tr(W(a)g(a)) = 0. 
For different a we generate randomly states g r with 
\\g(a) — g r \\ < 0.2 in Hilbert Schmidt norm ^(|. For the 
entangled states in this ball we determine the probability 
of detecting it via the witness W(a) and via a nonlinear 
criterion. First, we consider J 7 ^ as in Eq. ©, using \ip 2 ) 
and 2-^2. Then, we consider as in Eq. using 
all nonlinearities. We also test the criterion in Eq. J7J), 
corresponding to all possible T 1 ^ 1 . Finally, we investigate 
JF( 3 ) from Eq. ||ÜJ using ^3), |"04) and — c 2 . The re- 
sults are shown in Fig. 1. Clearly, the nonlinear criteria 
improve the witness significantly [17j . 

Let us now discuss generalizations of our approach 
to other entanglement witnesses, which are not related 
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FIG. 1: Probability of detecting a state with a witness and 
with nonlinear criteria. See text for details. 

to the PPT criterion. This can be done via positive 
maps [Üj . Let Hb and Hc be Hilbert spaces and let 
B{TLi) denote the linear operators on it. A linear map 
A : B{H B ) -> B(H C ) fulfiUing A(X) > for X > and 
A(X>) = A.(Xy is called positive (P). A positive map is 
completely positive (CP) when for an arbitrary TLa the 
map \a ® A is P, otherwise, it is positive, but not com- 
pletely positive (PnCP). Here, IL4 denotes the identity on 
B(H.a)- For example, the transposition is PnCP: While 
X > implies X T > the partial transposition does not 
preserve the positivity of a state. 

Indeed, it has been shown 0, 0] that a state g € 
B(TÍa) ® B(Hb) is separable iff for all P maps A the 
relation ILa ® A(g) > 0. holds. Consequently, if g is en- 
tangled there must be a trace decreasing PnCP map A 
where I J 4<g)A(^) has a negative eigenvalue A_ and a corre- 
sponding eigenvector \<f>). Taking (Ia<8>A) + as the adjoint 
of the map (ILa <g> A) with respect to the scalar product 
Tr(X^Y) a witness detecting g is given by 

W = (]U®A)+(|<£)<0|), (14) 

since we have Tr[pW) = Tr[p(I A O A) + (\4>) {(j>\)] = 
Tt[Ia ® A(g)\4>) = A_. Rcplacing the partial trans- 
position by the map (1a <8> A) + this witness can then be 
improved as in Observation 1. 

For an arbitrary witness, we make use of the 
Jamiolkowski isomorphism 0, between operators 
and maps. According to this, an operator E on B(Hb) <8> 
B(Hc) corresponds to a map e : B(Hb) — * B(TLc) act- 
ing as e{g) = TrB{Eg T ® 1c). Conversely, we have E — 
{lB>®e)(\è+)((j)+\) where Hb>=Hb and \<j>+) = Ei l**> 
is a maximally entangled state on H.B 1 <8 7~Íb- The key 
point is that if E is an entanglement witness, then e is a 
PnCP map [ï^. |2(ÏI|. Hence, any witness can be written 
as in Eq. (|14f> and finally we have: 

Theorem 1. Any bipartite entanglement witness can 
be improved by nonlinear corrections using the methods 
of Observation 1. 
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In conclusion, we have shown that all bipartite entan- 
glement witnesses can be improved by nonlinear expres- 
sions. These nonlinear witnesses are straightforward to 
calculate and can also be directly implcmentcd in exper- 
iments, often without any extra effort. It is tempting to 
extend these constructions to the multipartite scenario. 
Here, this challenge remains undone. 

We thank H.J. Briegel, M. Curty, J. Eisert, P. Hyllus, 
B. Kraus, M. Lewenstein, M. Piani and G. Tóth for dis- 
cussions. This work has been supported by the FWF, the 
DFG (Emmy-Noether-Programm) and the EU (OLAQI, 
PROSECCO, QAP, QUPRODIS, SCALA). 

Appendix — Here we prové Observation 2. (a) First, 
note that Y — Q TB \(f)){<j)\Q TB =■ \ri){f}\ is projector 
onto a not normalized state I77). To derive a crite- 
rion for the strength of J 7 ^ 1 ' we have to maximizc 
(A TB )((A TB )t)/ s (V>) = (ip\r))(r)\ip)/s(<ilj) for a given 
over all \ip). Lct \rj) = J2i a>i\H) and \tp) = Ei be the 
Schmidt decompositions of and \ip), with decreasingly 
ordered Schmidt coefficients. We first show that for fixed 
Schmidt coefficients bi it is optimal to take \ii) = \ii). 
We have |(^)[= |E« *M*IÍÏ>I = E,, 
where U^- = (i\j} A and = (i\j} B are unitary. Defin- 
ing üij = y/aïy/~bjU,{j and bij — \fai-\/bjU^ this can 
be written as a sealar product which is maximal, when 
a^, bij are paral·lel. Due to the Cauchy-Schwarz in- 
equality, we can then assume without loosing general- 
ity < Ey^jl^tíl 8 - Since U$ is unitary, \U$\ 2 

is doubly stochastic, i.e. its row and column sums 
equal one. Due to Birkhoff's theorem 21] it can be 
written as a convex combination of permutation matri- 
ces. For a permutation ir we have Ei a i^U) < Ei a i^i 
due to the ordering of the eti and bj. This implics fi- 
nally |(?7|V·')I ^ E,- a ibj where equality is achieved when 
= \ii). Having fixed — we only have max- 
imize Ei(^/^o)aj over all bi. This is maximal when 
(bi/bo) = 1 for all i, which corresponds to a maximally cn- 
tangled = l/v^Ej I**)- But then {ip\rj) (ri\ip) / s(ip) = 
(Ei a ï) 2 corresponds to the right hand side of Eq. J7J. 
(b) First, note that for an arbitrary basis \if>i) in Eq. © 

E í =i(^ Ts )((^ Ts ) t > = Tr{Y) = mg^m. This 
shows Eq. Jfjj. (c) Assume that g fulfills Eq. JSJ. Since 
Y is of rank one, there is a single \tp') such that Tr(Y) = 
(tp'\Y\ip'). So, if we take X = \<f>){íp'\ then the witness 
T = (\(/>)(<t>\ TB ) - (A T «)((A T «)t) deteets g. The wit- 
ness P = {\4>){4>\ Tb ) - {X Tb ){{X Tb )^)/s{^) is stronger, 
and deteets it as well. Finally, we have to show that 
Eq. |0 is never vàlid for a PPT state g. In view of 
the proof of (a), it suffices to show Q := (<fi\g TB \<j>) — 
(4>\g TB \<t>){<j>\g TB IV , )/ s ('0) > for arbitrary \<f>) and max- 
imally entangled l^). Since g Tu > we can define R = 

VOT^MVIF and S = VT^WWv^/sW, 
then Q = Tr(R) - Tr(RS) = Tr{R{t - S)). Since 
R > and S > it suffices to show that Tr(S) < 1, 
then (1 - S) > follows and finally Q > 0. We have 



Tr(S) = (ip\g TB / s(4>) ■ Now, we use the known fact 
that a witness like W = s(ip)l — \tp)(ip\ where \ip) is 
maximally entangled, deteets no PPT states 0. Since 
Q T B ig ppTj thig implies tliat Tr(S) < 1. □ 
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